## Loading required package: foreach
## 
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
## 
##     accumulate, when
## Loading required package: iterators
## 
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
## 
##     hour, isoweek, mday, minute, month, quarter, second, wday, week,
##     yday, year
## The following objects are masked from 'package:dplyr':
## 
##     between, first, last
## The following object is masked from 'package:purrr':
## 
##     transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.072208 -0.006215 -0.002112  0.005087  0.048581 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0101285  0.0002771   36.56   <2e-16 ***
## test_set$FCH4 0.2312130  0.0097955   23.60   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009789 on 1686 degrees of freedom
## Multiple R-squared:  0.2484, Adjusted R-squared:  0.2479 
## F-statistic: 557.1 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.168184 -0.004063 -0.000949  0.003350  0.090493 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063687  0.0001099   57.96   <2e-16 ***
## result$FCH4 0.5220989  0.0040188  129.91   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00779 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7143, Adjusted R-squared:  0.7143 
## F-statistic: 1.688e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060437 -0.006980 -0.001889  0.004907  0.052850 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0107836  0.0002939   36.69   <2e-16 ***
## test_set$FCH4 0.2188734  0.0105320   20.78   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01049 on 1686 degrees of freedom
## Multiple R-squared:  0.2039, Adjusted R-squared:  0.2035 
## F-statistic: 431.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.159933 -0.004229 -0.000997  0.002916  0.107940 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0066272  0.0001145    57.9   <2e-16 ***
## result$FCH4 0.5211818  0.0041861   124.5   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008114 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.6967, Adjusted R-squared:  0.6966 
## F-statistic: 1.55e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.049790 -0.006012 -0.001369  0.005024  0.066616 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0097941  0.0002695   36.34   <2e-16 ***
## test_set$FCH4 0.2580954  0.0105396   24.49   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009462 on 1686 degrees of freedom
## Multiple R-squared:  0.2624, Adjusted R-squared:  0.2619 
## F-statistic: 599.7 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.124486 -0.003897 -0.000916  0.003160  0.124550 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062906  0.0001051   59.83   <2e-16 ***
## result$FCH4 0.5428552  0.0038450  141.19   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007453 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7471, Adjusted R-squared:  0.747 
## F-statistic: 1.993e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.044136 -0.006188 -0.001447  0.004768  0.067381 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0097693  0.0002867   34.08   <2e-16 ***
## test_set$FCH4 0.2705066  0.0108225   25.00   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01004 on 1686 degrees of freedom
## Multiple R-squared:  0.2704, Adjusted R-squared:  0.2699 
## F-statistic: 624.7 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.078096 -0.004011 -0.000886  0.003129  0.115947 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062343  0.0001073   58.13   <2e-16 ***
## result$FCH4 0.5442627  0.0039225  138.75   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007603 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7404, Adjusted R-squared:  0.7404 
## F-statistic: 1.925e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060971 -0.006431 -0.001339  0.005305  0.060552 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0097295  0.0002955   32.92   <2e-16 ***
## test_set$FCH4 0.3161960  0.0117849   26.83   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0102 on 1686 degrees of freedom
## Multiple R-squared:  0.2992, Adjusted R-squared:  0.2988 
## F-statistic: 719.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.080506 -0.004035 -0.000927  0.003174  0.085621 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.006229   0.000104    59.9   <2e-16 ***
## result$FCH4 0.557341   0.003803   146.5   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007372 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7609, Adjusted R-squared:  0.7608 
## F-statistic: 2.147e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.042142 -0.007056 -0.002321  0.005192  0.090158 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0107851  0.0003239   33.30   <2e-16 ***
## test_set$FCH4 0.2559779  0.0114887   22.28   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01163 on 1686 degrees of freedom
## Multiple R-squared:  0.2275, Adjusted R-squared:  0.227 
## F-statistic: 496.4 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.097414 -0.004061 -0.001158  0.002898  0.132175 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062451  0.0001204   51.87   <2e-16 ***
## result$FCH4 0.5599477  0.0044035  127.16   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008536 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7055, Adjusted R-squared:  0.7055 
## F-statistic: 1.617e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.102662 -0.006545 -0.001355  0.005129  0.046061 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0105528  0.0002852    37.0   <2e-16 ***
## test_set$FCH4 0.2198590  0.0102740    21.4   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01015 on 1686 degrees of freedom
## Multiple R-squared:  0.2136, Adjusted R-squared:  0.2131 
## F-statistic: 457.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.225351 -0.004068 -0.000969  0.003185  0.079733 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0064394  0.0001128   57.11   <2e-16 ***
## result$FCH4 0.5270819  0.0041237  127.82   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007993 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7077, Adjusted R-squared:  0.7076 
## F-statistic: 1.634e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.071279 -0.006088 -0.001542  0.004459  0.051965 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0101719  0.0002817   36.11   <2e-16 ***
## test_set$FCH4 0.2509637  0.0105253   23.84   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00981 on 1686 degrees of freedom
## Multiple R-squared:  0.2522, Adjusted R-squared:  0.2517 
## F-statistic: 568.5 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.164506 -0.003931 -0.000896  0.003057  0.091802 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063248  0.0001078   58.67   <2e-16 ***
## result$FCH4 0.5340847  0.0039424  135.47   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007642 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7311, Adjusted R-squared:  0.7311 
## F-statistic: 1.835e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.053052 -0.006865 -0.002164  0.004994  0.062974 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0103832  0.0002952   35.17   <2e-16 ***
## test_set$FCH4 0.2911413  0.0115847   25.13   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01029 on 1686 degrees of freedom
## Multiple R-squared:  0.2725, Adjusted R-squared:  0.2721 
## F-statistic: 631.6 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.062335 -0.004168 -0.001128  0.003099  0.106795 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063393  0.0001069    59.3   <2e-16 ***
## result$FCH4 0.5563607  0.0039096   142.3   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007578 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:   0.75,  Adjusted R-squared:   0.75 
## F-statistic: 2.025e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.086313 -0.007140 -0.001803  0.005344  0.053994 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0113257  0.0002863   39.56   <2e-16 ***
## test_set$FCH4 0.1865034  0.0096318   19.36   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01038 on 1686 degrees of freedom
## Multiple R-squared:  0.1819, Adjusted R-squared:  0.1814 
## F-statistic: 374.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.211087 -0.004384 -0.001055  0.003223  0.115594 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0068941  0.0001183   58.27   <2e-16 ***
## result$FCH4 0.4995486  0.0043271  115.45   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008387 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.6638, Adjusted R-squared:  0.6638 
## F-statistic: 1.333e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.061838 -0.005811 -0.001497  0.004491  0.058570 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0099197  0.0002776   35.74   <2e-16 ***
## test_set$FCH4 0.2550037  0.0107117   23.81   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009777 on 1686 degrees of freedom
## Multiple R-squared:  0.2516, Adjusted R-squared:  0.2511 
## F-statistic: 566.7 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.136980 -0.004027 -0.000989  0.003057  0.091629 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062858  0.0001073   58.61   <2e-16 ***
## result$FCH4 0.5448636  0.0039225  138.91   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007603 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7409, Adjusted R-squared:  0.7408 
## F-statistic: 1.929e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.048217 -0.006349 -0.001687  0.005078  0.054675 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0100743  0.0002925   34.44   <2e-16 ***
## test_set$FCH4 0.2773132  0.0115507   24.01   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01009 on 1686 degrees of freedom
## Multiple R-squared:  0.2548, Adjusted R-squared:  0.2543 
## F-statistic: 576.4 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.085364 -0.003945 -0.000902  0.003066  0.093297 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0061941  0.0001067   58.07   <2e-16 ***
## result$FCH4 0.5533471  0.0039011  141.84   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007562 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7488, Adjusted R-squared:  0.7488 
## F-statistic: 2.012e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.085089 -0.006953 -0.001988  0.004998  0.071706 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0099375  0.0003258    30.5   <2e-16 ***
## test_set$FCH4 0.2958720  0.0124341    23.8   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01123 on 1686 degrees of freedom
## Multiple R-squared:  0.2514, Adjusted R-squared:  0.251 
## F-statistic: 566.2 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.180234 -0.004005 -0.001030  0.003031  0.119398 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0057118  0.0001143   49.95   <2e-16 ***
## result$FCH4 0.5856925  0.0041817  140.06   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008106 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.744,  Adjusted R-squared:  0.744 
## F-statistic: 1.962e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060942 -0.006939 -0.001814  0.005275  0.074804 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.010712   0.000304   35.23   <2e-16 ***
## test_set$FCH4 0.233880   0.011623   20.12   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01074 on 1686 degrees of freedom
## Multiple R-squared:  0.1937, Adjusted R-squared:  0.1932 
## F-statistic: 404.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.133127 -0.004161 -0.000960  0.003293  0.122737 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0064488  0.0001131   56.99   <2e-16 ***
## result$FCH4 0.5418519  0.0041381  130.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008021 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7176, Adjusted R-squared:  0.7175 
## F-statistic: 1.715e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.067802 -0.006469 -0.001742  0.004936  0.063292 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0106747  0.0002769   38.55   <2e-16 ***
## test_set$FCH4 0.2352920  0.0101966   23.08   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009838 on 1686 degrees of freedom
## Multiple R-squared:   0.24,  Adjusted R-squared:  0.2396 
## F-statistic: 532.5 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.164741 -0.004099 -0.000996  0.003139  0.123711 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0065806  0.0001095   60.08   <2e-16 ***
## result$FCH4 0.5299611  0.0040058  132.30   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007765 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7217, Adjusted R-squared:  0.7217 
## F-statistic: 1.75e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.089485 -0.006665 -0.001608  0.005269  0.056685 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.010933   0.000285   38.37   <2e-16 ***
## test_set$FCH4 0.199724   0.009266   21.55   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01038 on 1686 degrees of freedom
## Multiple R-squared:  0.216,  Adjusted R-squared:  0.2156 
## F-statistic: 464.6 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.205545 -0.004388 -0.001019  0.003426  0.104318 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0069286  0.0001172    59.1   <2e-16 ***
## result$FCH4 0.4906196  0.0042878   114.4   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008311 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.6599, Adjusted R-squared:  0.6598 
## F-statistic: 1.309e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.061226 -0.006867 -0.001887  0.005102  0.056247 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0105388  0.0003071   34.32   <2e-16 ***
## test_set$FCH4 0.2756657  0.0118356   23.29   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01064 on 1686 degrees of freedom
## Multiple R-squared:  0.2434, Adjusted R-squared:  0.243 
## F-statistic: 542.5 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.084846 -0.004159 -0.001164  0.003180  0.091598 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063361  0.0001097   57.75   <2e-16 ***
## result$FCH4 0.5490025  0.0040127  136.82   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007778 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.735,  Adjusted R-squared:  0.735 
## F-statistic: 1.872e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.113537 -0.006393 -0.001591  0.005115  0.057936 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0106040  0.0002985   35.52   <2e-16 ***
## test_set$FCH4 0.2354589  0.0111378   21.14   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01047 on 1686 degrees of freedom
## Multiple R-squared:  0.2095, Adjusted R-squared:  0.2091 
## F-statistic: 446.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.232980 -0.004157 -0.000997  0.003369  0.072216 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0064144  0.0001125   57.01   <2e-16 ***
## result$FCH4 0.5349993  0.0041151  130.01   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007976 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.7147, Adjusted R-squared:  0.7146 
## F-statistic: 1.69e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.073614 -0.006999 -0.001671  0.005542  0.066899 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.010556   0.000313   33.72   <2e-16 ***
## test_set$FCH4 0.254681   0.010930   23.30   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01111 on 1686 degrees of freedom
## Multiple R-squared:  0.2436, Adjusted R-squared:  0.2431 
## F-statistic: 542.9 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.173372 -0.004090 -0.000957  0.003194  0.116958 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0060873  0.0001174   51.85   <2e-16 ***
## result$FCH4 0.5586618  0.0042938  130.11   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008323 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.715,  Adjusted R-squared:  0.7149 
## F-statistic: 1.693e+04 on 1 and 6749 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.103887 -0.006922 -0.002189  0.005336  0.060620 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0109562  0.0002942   37.24   <2e-16 ***
## test_set$FCH4 0.2187464  0.0100288   21.81   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01048 on 1686 degrees of freedom
## Multiple R-squared:  0.2201, Adjusted R-squared:  0.2196 
## F-statistic: 475.8 on 1 and 1686 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.221616 -0.004182 -0.001106  0.003307  0.088614 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0065925  0.0001145   57.56   <2e-16 ***
## result$FCH4 0.5145586  0.0041891  122.83   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00812 on 6749 degrees of freedom
##   (7936 observations deleted due to missingness)
## Multiple R-squared:  0.6909, Adjusted R-squared:  0.6909 
## F-statistic: 1.509e+04 on 1 and 6749 DF,  p-value: < 2.2e-16